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The canonical quantum Hamiltonian eigenvalue problem for an anhar- 
monic oscillator with a Lagrangian L = 0^/2 — m?(fp' /2 — gm^cj)'^ is nu- 
merically solved in two ways. One of the ways uses a plain cutoff on the 
number of basis states and the other employs a renormalization group pro- 
cedure. The latter yields superior results to the former because it allows 
one to calculate the effective Hamiltonians. Matrices of effective Hamilto- 
nians are quite small in comparison to the initial cutoff but nevertheless 
yield accurate eigenvalues thanks to the fact that just eight of their highest- 
energy matrix elements are proper functions of the small effective cutoff. 
We explain how these cutoff-dependent matrix elements emerge from the 
structure of the Hamiltonian and the renormalization group recursion, and 
we show that such small number of cutoff-dependent terms is sufficient to 
renormalize any band-diagonal Hamiltonian. 



PACS numbers: ll.lO.Gh 



1. Introduction 



The purpose of this article is to show that a renormalization group 
procedure (RGP) of Hamiltonians, based on Gaussian elimination of ba- 
sis states [1], provides a convenient approach to the Hamiltonian matrix 
diagonalization problem for an elementary oscillator with quartic anhar- 
monicity, 

= 1 + !I^ + ,,„V, (1) 

where 5 > is a dimensionless coupling constant and the units are chosen 
such that h = 1 and c = 1. By "elementary" we mean that we consider (j) 
a single variable, describing only one mode of a scalar field. Thus, instead 
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of a quantum field theory problem, we have a simple quantum mechanics 
problem. The utility of RGP results from the band-diagonal structure of 
the oscillator's Hamiltonian. Having written the initial Hamiltonian in the 
second-quantized form, 



(the irrelevant constant +m/2 was suppressed), the Hamiltonian matrix is 
obtained by evaluating matrix elements in the basis of normalized eigen- 
states of Hq = ma) a. 

In general, band-diagonal structure of a Hamiltonian corresponds to the 
situation, when an interaction cannot mix the states belonging to distinct 
energy scales. Such an interaction is much easier to understand then the one 
which does not possess this property. For this reason, building a unitary 
transformation which decouples states distinct in energy from each other 
may be the main part of solving a Hamiltonian eigenvalue problem involv- 
ing interactions of the latter kind. The examples of methods dealing with 
such unitary decoupling are the similarity renormalization group (SRG) 
scheme O [3] and (closely related) Wegner's flow equation [HEIE]. Both 
methods bring (arbitrarily chosen) initial Hamiltonian to a band-diagonal 
form. The oscillator Hamiltonian possesses a band diagonal structure from 
the beginning, i.e., without any need for using SRG to make it so. Thus, 
this paper can be considered as dealing with the second stage of solving the 
eigenvalue problem after the band diagonal structure is already achieved. 

Namely, in this paper we address the problem of precise and systematic 
numerical calculation of small eigenvalues of the band-diagonal Hamiltonian, 
using the one in Eq. ([2]) as an example. We show that the Wilsonian RGP in 
the case of any band-diagonal Hamiltonian is reduced to a recursive relation 
for only a few (six in the considered example) independent Hamiltonian 
matrix elements. Thus, what we do here may be considered a convenient tool 
for analysing band-diagonal Hamiltonians emerging from aforementioned 
methods. Further discussion of this possibility in the case of realistic field- 
theoretic Hamiltonians can be found in [3], see Appendix A there. 

In zero-dimensional oscillator there are no ultraviolet divergences, so 
the renormalization issues in this model are limited to finite cutoff depen- 
dence. Such finite cutoff dependence strictly disappears from the smallest 
eigenvalues when the cutoff is sent to infinity. Thus, there is no need to 
introduce diverging counterterms. Nevertheless, RGP is still of great value, 
because it gets rid of the spurious finite cutoff dependence of small Hamil- 
tonian eigenvalues, even for very small cutoffs. The numerical recipe is very 
fast, precise and nonperturbative. So, it is valid for all values of g. The 
low-energy spectra of the Hamiltonian ^ obtained with RGP method will 
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be compared with the corresponding spectra obtained by simply cutting off 
the initial Hamiltonian at the same small cutoff, for different values of g. 

RGPs are extensively used as numerical tools since Wilson has delivered 
the algorithm for treating the systems in which Kondo effect |8j appears. 
The algorithm considered in this paper resembles the one considered in [7], 
but does not include a logarithmic discretization of any conduction band. It 
does, however, deal with a discrete set of states that is similar to an infinite 
set of intervals of pion momentum in [9] and a discrete set in [1] . Ref . |10] 
provides an overview of applications of numerical RGPs to the quantum 
impurity systems. Here we show that an appropriate numerical RGP can 
be always applied to any quantum Hamiltonian that is band-diagonal. In 
particular, we have in mind systems of particles for which SRG renders a 
Hamiltonian of the band-diagonal form [3] . 

Our approach is not the only one possible in the case of a quartic oscilla- 
tor. Another approach, based directly on SRG, was presented in [11]. In the 
case of oscillator, the procedure of Ref. [11] is simple and an effective gen- 
erator of SRG transformation is found using specific simplifications. Most 
of the work is done analytically. However, for an arbitrary band-diagonal 
Hamiltonian such simple generator does not exist and numerical calcula- 
tions based on SRG are complex. In fact, in such circumstances they are 
less effective then these based on Wilsonian RGP that are presented in this 
article using the example of quartic oscillator. 

The paper is organized as follows. The RGP carried out in this article 
is described in Section [21 The resulting evolution of Hamiltonian matrix el- 
ements is described in Section [3j Section [5] compares results obtained using 
RGP with results obtained using a plain cutoff (PC) of the Hamiltonian ma- 
trix, without inclusion of any corresponding changes in its matrix elements. 
The PC procedure has no a priori justification for general Hamiltonians but 
works for large enough cutoffs for band-diagonal Hamiltonians, see Section [5] 
for details. The effective Hamiltonians we obtain are described in Section [5j 
Section [6] concludes the article with a summary of accuracy achieved using 
the RGP. 

2. Description of RGP carried out in this article 

We set m = 1 and denote Hj = [H — Hq)/ g. In the basis of normalized 
eigenstates of Hq^ \k) = (A;!)~-^/^(a^)''|0), H has matrix elements 



Hki = {k\H\l) 



= [k + ^g{2k^ + 2k + l)]5ki + 

+ g{Al + Q)^{l + l){l + 2)5kii+2) + 
+ g{ik + 6)^/{kTr){k + 2)5i^k+2) + 
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+ gVik + l){k + 2){k + ?,){k + 4)Ji(fc+4) . 

Due to the Kronecker delta functions most of the matrix elements are zero 
and the Hamiltonian matrix is band-diagonal with only five near-diagonals 
different from zero. 

As the first step of RGP, we define the Hamiltonian with a big cut- 
off i.e., we reduce the space of states to the one spanned by the set 
{|0), |1), . . . , |iV)} and write the Hamiltonian as (A^ + 1) x (A^ + 1) matrix 
with elements given by Eq. ([3]). 

Then, we denote {k\il)) = tpk and write the resulting Hamiltonian matrix 
eigenvalue equation in the form 

N 

Y^Hki^i = Ei^k- (4) 

1=0 

We extract 'i/'TV from the highest energy equation, A; = A^, in the set ^ of 
A^ + 1 equations, and use it in the remaining A^ equations in the set @ with 
k < N. This is essentially Gaussian elimination. The resulting set of only 
A^ equations is 

Taking into account that H = Hq + gHj, we arrive at the equation 

j:(Ho,i+gH,,, + g \ = ^^^^ (6) 

T~i \ - \n-ONN + gtll NN) J 
L-0 V ^ . 

which is an eigenvalue problem for a new matrix with elements i^ofei + 
gHj^^i^{E). We call it the matrix of effective Hamiltonian. It can be used 
as a starting point for the next Gaussian elimination step, but we notice 
that there appears a difficulty: Hj^^ depends on an unknown eigenvalue E. 
However, we are interested only in small eigenvalues of H, so small that 
they are negligible when compared with Hq^n + gHj nj^. The latter is the 
highest-index diagonal matrix element of H that appears in the denomi- 
nator on the left-hand side in Eq. ([6]). Thus, we bypass the difficulty by 
introducing the approximation 

E — (Hqnn + gHj nn) ~ —{HoNN+gHiNN)- (7) 
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This approximation limits our ability to precisely calculate eigenvalues E 
using our RGP to lowest eigenvalues. For this price we can do as many 
Gaussian steps as we like, using the relation 

Hiki -Hj,i-g- (,) , (8) 

Honn + g^Inn 

with n = N—j, provided that we keep satisfying condition ([7]) for every step 

number j. We may suppose that for some M = N — j, Hqmm + 9Hj^j^^ « 
E and the condition ([7|) may be eventually violated. On the other hand, 
we can continue a calculation assuming the simplification and using the 
approximate Gaussian steps to eliminate states |A^), \N — 1), . . . , |n + 1) 
down to a small n. We analyse results of such steps for a set of smallest 
eigenvalues for which the violation of our approximation appears quite small. 

From now on, the Hamiltonian reduced in size by performing N — n 
approximate Gaussian steps (in other words the renormalized Hamiltonian) 
is denoted by H^'^{n). On the other hand, a Hamiltonian matrix reduced 
to the same size with a PC is denoted by Hf^{n). This means that both 
Hf'^in) and Hf^{n) are (n+1) x (n+1) matrices but their matrix elements 
are different. The matrix elements of H^^{n) are to be calculated and the 
matrix elements of Hf'~"{n) are directly given by the right-hand side of 
Eq. Q as the terms proportional to g. 

Note that some matrix elements of the renormalized matrix depend on 
N and n. If the dependence on N were divergent or otherwise significant, 
it would have to be removed by counterterms in H^'^\ Since in our model 
there are no divergences, and the limit of — >■ oo is easily achieved for small 
eigenvalues, one can simply use some big N and calculate the corresponding 
matrices of operators H^^{n). N sufficiently large for working without 
counterterms to be valid, is determined in Sec. HI H^^{n) will have some 
matrix elements that vary with n. The only variation we study here is the 
one obtained assuming condition ([7|). A more precise study than the one 
described here would be required to identify consequences of a finite ratio 

E/Hnn- 

3. RGP evolution of Hamiltonian matrix elements 

3.1. General hand- diagonal Hamiltonian 

Matrix elements of interaction part of any band-diagonal Hamiltonian 
can be written in the form 

m 

Hiki = hi{k)5k-i,i (9) 
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for some integer m. From Eq. ([8]) follows that H^SiN — 1) / Hf^i{N) = 
Hf^i{N) only for such k that Hf^^{N) / and Hfg{N) ^ 0. This implies 
both k,l > N — m. Thus, after first Gaussian step, done in appoximation 
([7]) or exactly, H^'-' is still band-diagonal. Moreover, only matrix elements 
situated in m x m submatrix in highest-energy corner of H^^{N — 1) can 
be different then corresponding matrix elements of Hf'-'{N — 1). We see, 
that we can repeat this reasoning recursively, with N changed into n (the 
highest index after some number of Gaussian steps). 

Thus, we can conclude that H^^{n) is band-diagonal for any n, and that 
not more then its matrix elements, located in the high-energy corner, 
does really depend on n. Eq. ([8]) allows to write down the set of (conjugated) 
recursions for these matrix elements. Numerical iteration of such a set can 
be easily done using computer. 



3.2. Oscillator example 

First, we see that in the case of oscillator's Hamiltonian ([3]), m = 4. 
Further, from Eqs. ^ and ([8|) we see, that no new nonzero matrix elements 
appear in the 4x4 submatrix in the high-energy corner during the Gaussian 
step. Thus, only eight matrix elements depend on n. Finally, the symmetry 
of Hamiltonian matrix is preserved by RGP transformation, so there are only 
six independent functions of n in H^^{n). To number them, we introduce 
a function 



i{k,l) 



and denote 



j for 

5 for 

6 for 



k = l = n- {j-l) and j e {1,2,3,4}, 
{k,l) = (n,n — 2) or {k,l) = (n — 2,n), 
{k,l) = (n — l,n — 3) or {k,l) = {n — 2>,n- 



Hfi{n) 



(10) 



(11) 



This mean that ^j(n=A^) = 1 and the deviation of Ci{k,i)i''T') from unity is the 
relative correction to Hf^^ resulting from RGP. Values of g^i are effective 
coupling constants in the highest-energy 4x4 submatrix of H^^{n). 

From Eq. ([HD, we obtain a set of recursions for coefficients (,i{n). In the 
case of oscillator six conjugated first order recursive equations can be simply 
split into two sets of three equations of the second order. One such set is 



Un-2) = 6(n)-/°2 



e3(n-2) 



1-5 



/22 n + gfooCi{n) - E' 

/o^4 1 

/44 n + 5/006 {n) - E' 



(12) 
(13) 
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i^{n-2) = l-g— ■ — „ . . . ^, 14 

/24 n + 5/006 (f^j - ^ 

where fki = are functions of n, explicitly known from Eq. 

The second set is almost identical. The only differences are that all down 
indices are increased by one and (n — 1) appears instead of n in the de- 
numerator on the right-hand side. For large n the two sets evolve nearly 
identically. It is evident that the coefficients ^j(n) are easy to calculate 
numerically. 



4. Comparison of RGP and PC 

In Eq. (I6|), we can see that Hf^{n) can be interpreted as Hf^{N) after 
N — n Gaussian steps performed in the approximation 

Hjki + g— w Hjki- (15) 

-C/ — rlNN 

This approximation neglects terms obtained in the approximation ([7|) . There- 
fore, we can predict, that the results obtained using the RGP will be more 
accurate then the ones obtained using the PC. 

It is not clear how to compare the accuracies of the two procedures, since 
we don't know the exact eigenvalues of the infinite Hamiltonian. However, 
in our simple model, we can very precisely calculate small eigenvalues of the 
Hamiltonian numerically diagonalizing H^^{M) for some big cutoff M. 

In all computations, the numbers were handled with precision of 53 bi- 
nary places (almost 16 decimal places). The numerical calculations show 
that for M G [200, 1000] and g G {0.01, 0.1, 1, 10} the three smallest eigen- 
values, denoted Eq, Ei, E2, do not depend on M with a relative precision 
10"^°. For this reason, we can consider them to be accurate results and 
define the (relative) accuracy of RGP and PC as \E(^^{n) - Ei\/Ei and 
|£'f'^(n) — Ei\/Ei, respectively. Also for this reason N = 200 will be con- 
sidered sufficiently large initial cutoff and it will be used as a starting point 
for RGP described in previous sections. 

We shall illustrate our results on examples with three different values 
of the initial coupling constant: g = 0.01, g = 1 and g = 10. These 
values are chosen because each one illustrates a different regime of RGP 
results, g = 0.01 illustrates how RGP works in the case of interactions 
which only slightly change eigenvalues of Hq. (7 = 10 illustrates how RGP 
works in the case of interactions dominating the system (which manifest 
itself in eigenvalues much different then those of Hq). g = I illustrates the 
intermediate situation. The results are discussed separately for different 
values of g. 
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4.1. Results for g = 0.01 

For g = 0.01, the accurate eigenvalues are Eq = 0.1687726041, Ei = 
1.716925770 and E2 = 3.602838696. Both procedures, RGP and PC, give 
precise results. Even for n as small as n = 10, Eq is recovered with accuracy 
4 • 10-1° and Ei and E2 are recovered with accuracy 2-10 ^ in the case of 
RGP. In the case of PC, the accuracy of these three eigenvalues is about 
10""^. We see that RGP increases accuracy of predicting £^0 about 500 times. 
Also E^'^ and are calculated about 5 times more precisely then E^'-' 
and E2^ . Nevertheless, accuracy obtained with PC is also brilliant. 

4-2. Results for g = 1 

For g = 1 and n = 10, eigenvalues Ef-" and i?!^*^ differ from eigenvalues 
E^'^ and £^^*^ much more then for g = 0.01. We can say that the difference 
is qualitative. The correct eigenvalues are Ei = 3.521565666 and E2 = 
7.263980184. E^'^ has accuracy 3%, while E^'^ only 21%. For Ei, the 
difference is 0.25% to 6% in favor of RGP. 

The correct ground state energy is £^0 = 0.6487889141. The errors are 
2-10"^ in RGP and 6-10-^ in PC. Both values can be considered satisfactory 
if the required accuracy is of the order of 1%. 

4.3. Results for g = 10 

Utility of RGP is most clearly visible in the case of g = 10. This is giant, 
but still realistic coupling. For example the pion-nucleon coupling constant 
is about 13 [12]. Moreover, there is no a priori reason for values of couplings 
emerging from SRG procedure to be small. 

The accurate eigenvalues in this case are Eq = 1.826275924, Ei = 
7.790412053 and E3 = 15.70695963. The results for eigenvalue Eq, obtained 
with RGP and with PC, as a function of cutoff n are plotted in Fig. [TJ One 
sees that PC with highest index smaller than ~15 produces a huge error in 
the eigenvalue. In the same cutoff region, RGP still gives a reasonable value 
oIEq. 

For n = 10, Eq^ has accuracy 0.35%, which is about 100 times better 
than 33% obtained with PC. These results prove that as long as E <^ Hnn ~ 
v?, RGP gives a reasonable estimation of an eigenvalue. Even for n = 4, 
RGP reproduces Eq with 5.5% accuracy while PC introduces about 230% 
of error. 

For n = A, the eigenvalues Ei and E2 are obtained with poor accuracy. 
In RGP, Ei^ and E2^ have errors 11% and 43%, respectively. Obtaining 
E^'-" and -E'^^'^ does not make any sense. 
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Figure 1. The ground state energy of H for g = 10 determined with RGP and PC, 
compared with the value achieved in the accurate numerical calculation, Num. 

5. Effective Hamiltonians 

The important question is how it happens that the RGP produces a 
better result then PC. In PC, there appears a non-physical parameter n. 
If n is not large enough, even small eigenvalues depend on n. The key 
idea behind the RGP is to base physical predictions (in this case the small 
eigenvalues of H) on the effective theory whose Hamiltonian Hf'~^{n) is 
calculated from the initial one, instead of plainly cutting off the initial theory 
to the arbitrary number of states. 

As it was explained in Sec. El in the case of band-diagonal Hamiltonian 
most of the matrix elements of the effective Hamiltonian stay invariant under 
the RGP transformation given by Eq. ([8|). 

The ratios of the changed terms in Hf^^{n) to the corresponding terms 
in Hf-" for g = 10 are presented in Fig. [21 One can see that for n < N 

< 6 ~ 6 < ~ < 6 ~ ^4 < 1 . (16) 

This means that these renormalized matrix elements are smaller then the 
original ones, but the sign remains unchanged. The approximate equalities 
correspond to the observation made in Sec. [3] that the six evolving matrix 
elements form two similar independent subsets. 

Within the 4x4 highest-energy corner of the Hamiltonian matrix, the 
closer is the matrix element to the cutoff corner, the bigger is the necessary 
correction. It is worth noticing that although considerably depart from 
unity, they still very weakly depend on n. Calculations for various g show 
that they also weakly depend on g. 
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Figure 2. The ratios of renormalized matrix elements to the original ones, as 
functions of n, obtained for the initial cutoff N = 200 and g = 10, see Eq. ([TT|). 



We have considered a simple version of RGP in application to an el- 
ementary Hamiltonian of a single-mode scalar field with an anharmonic 
interaction term proportional to ip^. We have calculated the ground state 
and two least excited states energies and compared results obtained using 
RGP to the results obtained using PC. 

Our results show that when the coupling constant g is large, the RGP is 
much more accurate than the PC procedure. For g = 10, RGP reproduces 
Eq with 5.5% accuracy in a 5 x 5 matrix, which can be considered a rea- 
sonable estimate. Such small matrices provide a model of small spaces of 
states that one can handle non-perturbatively using computers in realistic 
theories. In the matrix of the same size, PC gives the result with an enor- 
mous error (about 230%). RGP is significantly more accurate also for small 
values of g. 

In the presented oscillator model, the improvement due to RGP is 
achieved through a calculation of eight cutoff-dependent terms (only six 
of them being different) in the effective Hamiltonians with small cutoffs. A 
small number of terms is an unavoidable consequence of the band-diagonal 
structure of the initial Hamiltonian. This feature make RGP a poten- 
tially convenient tool for analysing band-diagonal Hamiltonians, in par- 
ticular these emerging form applying SRG to some realistic field-theoretic 
Hamiltonians. 



6. Conclusion 



The author would like to thank Stanislaw Glazek for many discussions. 
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